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Abstract 

■ The nuclear collective response is investigated in the framework 
, of a doorway picture in which the spreading width of the collective 
I motion is described as a coupling to more and more complex config- 

■ urations. It is shown that this coupling induces fluctuations of the 
, observed strength. In the case of a hierarchy of overlapping decay 
I channels, we observe Ericson fluctuations at different scales. Meth- 
ods for extracting these scales and the related lifetimes are discussed. 
Finally, we show that the coupling of different states at one level of 

Q I complexity to some common decay channels at the next level, may 

produce interference-like patterns in the nuclear response. This quan- 
tum effect leads to a new type of fluctuations with a typical width 
related to the level spacing. 



X 



1 Introduction 



Lifetimes and damping mechanisms are general questions in physics. In quan- 
tum mechanics, lifetimes are directly related to the width of the considered 
states through the Heisenberg relation. Lifetimes of simple systems, such 
as single particle states are often related to quantum tunnelling. In many- 
body systems, the damping of excitations appears more complex since it 
can involve many different processes. This is in particular the case for col- 
lective excitations of nuclei for which several origins for the observed width 
have been discussed so far[jl], ^ 0]. On the first hand, a fragmentation of the 
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collective response into several collective states is expected, leading to the so- 
called Landau damping. The collective motion can also be directly coupled 
to the continuum of escaping states[^ giving rise to the escape width. 
Finally, the collective strength can decay toward the compound nucleus con- 
figurations. This spreading width is often viewed in a doorway picture as 
the coupling to more and more complex states. For instance, a collective 
excitation built from particle-hole excitations can be coupled to two-particle 
two- hole states through the residual two-body interaction 0. Again, those 
states might themselves decay toward three-particle three-hole states because 
of two-body collisions. After many such steps up in complexity, this process 
eventually ends in the chaos of compound nucleus states[0, |], |10|( which 
may, in general, induce a fine-structure in the response of the nucleus |TT|]). 
In such a picture, collective excitations of many-body systems exhibit a large 
variety of time scales for the decay mechanism going from the short lifetime 
of the collective motion associated to a width of several MeV, to the long 
lived compound nucleus states with a typical width of the order of few eV or 
even less. 

Correlations are a well known tool to measure lifetime of quantal system. 
Indeed, when the level spacing is small compared with the typical width 
of each individual level the observed spectrum is known to exhibit Ericson 
fluctuations [1121. In such excitation spectra presents fluctuations with 



a correlation width characteristic of the average width of the overlapping 
states. As far as the collective response is concerned, because of the various 
scales involved in the decay process, one may wonder about the characteris- 
tics of the expected fluctuation pattern. Since the compound nucleus states 
can be considered as the true eigenstates of the many-body system, one 
would predict that the observed fluctuations should be characteristic of the 
associated lifetime. However, these fluctuations are over a so small scale that 
most experiments are unable to detect them. On the other hand, collective 
states have a large width and indeed poor resolution experiments exhibit 
resonant structure with several MeV width. Improving upon the resolution 
one thus expects to uncover more and more detailed structures directly re- 
lated to the different levels of complexity in the damping mechanism. The 
Landau spreading can be flrst studied. Then, according to the commonly 
admitted picture, one may look at the coupling to the two-particle two-hole 
states using a flner coarse graining, and so on and so forth down to the com- 
pound nucleus scale. Is it possible to observe such a multiscale structure in 
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the fluctuation pattern, from the collective mode toward different levels of 
complexity down to the compound nucleus chaos? In such a context, the 
response function may present more and more detailed structure when the 
resolution is improved, so one may even ask himself about a possible fractal 
nature of the response function. In this case, are the fluctuations self-similar 
or not? And so which interpretation can be given to the observed fluctuations 
function of the resolution in used? 
However, the above discussion might well be too simple. Indeed, the 
robustness of the various complexity levels against the effects of the residual 
interaction, has not been investigated. Even if from the classical point of 
view a typical path toward chaos may go through a sequence of bifurcations, 
there is no guarantee that, in quantum mechanics, the coupling of a regular 
collective motion to an ensemble of chaotic states follows this path. The 
simple picture of a hierarchy of more and more complex doorway states is 
typical of a perturbative approach, but in such large matrix, one strong 
element (compared with the typical level spacing) may be enough to change 
the spectrum and the properties of many eigenstates: i.e., at some point, the 
perturbative approach may break down. Moreover, many quantum effects 
can also be expected such as the possibility of interferences [O, O, [15 



or 



the indirect coupling of doorway states through their decay at higher orders 
in perturbative expansions. 

In this article we would like to present a critical discussion of the various 
questions raised above. The possibility of having multiscale fluctuations is 
flrst illustrated in section II. Since, it might be difficult to disentangle various 
scales in a fluctuating spectrum, we present in section III a brief discussion 
of the various methods for extracting these scales. Finally, in section IV, we 
demonstrate through analytical derivations and numerical calculations that 
fluctuations may also occur due to quantum interferences. 



2 Formalism and results 

Let us flrst briefly recall some known results about damping mechanism of 
doorway states and about the fluctuations of cross sections. In such a way, 
we will make clear the concepts and notations used in the article. 
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2.1 Ericson Fluctuations in doorway processes 



The collective response can be tested by applying an external field to the 
nucleus. If we assume that the ground state |0) is excited through a time- 
dependent field \D (^e^"^^^ + e^*^* j , where D is a hermitian operator. The 
linear response theory tells us that the state of the system isF][|T^ 

(1) 

where we have assumed that each eigenstate of the system has a finite 
lifetime^ = l/F^. Using the doorway notations D |0) = \CoU) and per- 
forming a Fourier transform of the observation D (t)), we get the 
corresponding spectral response which reads for positive energies 

^\{E-E,)+tT,/2 {E + E,) + iV,/2) 

Since the first term is dominant we can write R{E) ~ ^ \Coll)\^ / {E — E^ + iV ^/2). 
The associated strength function is given by: 

^ ' ^ ^ ^ " ^{E-E^f + rj/4 ^ ' 

When the various states are well separated, i.e. when the average 
distance AE" between two states is much larger than the individual width 
r^, the strength 5* [E) presents isolated peaks with a typical width V^. 
When the states are strongly overlapping, one can still extract information 
about the individual widths by looking at the fluctuations known as Ericson 



fluctuations [ll2|. Indeed, the overlap 0^= \Coll)\^ can be separated into 
two parts 

• an averaged part O {E^) which is a smooth function of E^. In the case 
of a resonance, this averaged part often take the shape of a gaussian or 



^in all the paper we used h = \ 

^ In the following, the states > are considered as eigenstates of an effective hamil- 
tonian. In the first part of this article the coupling of these degrees of freedom to the rest 
of the system will be simply taken into account as a finite lifetime F^. This is valid only 
for nonoverlapping decay channels as discussed in refs. |jl^, and in the last chapter. 
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a lorentzian centered on the collective energy with a width interpreted 
as the width of the doorway state [jl^ ; 



• a fluctuating one 50^ = 0^ — {E^). 

Then the response can also be splitted into two components (see ap- 
pendix 1) . On the one hand, it contains an averaged strength function 
S (E) proportional to O and on the other hand remains a fluctuating term 
6S (E) = S (E) — S (E) associated with 50^. These fluctuations still con- 
tain information about the averaged width T as it was demonstrated by 



Ericson[12 . To measure fluctuations, we can deflne the autocorrelation func- 



tion of the strength distribution as 

1 rEo+SE 

C{E) = — dE' SS {E') SS {E' + E) (4) 

oE JEo 

where SE is the averaging interval. If we assume random correlations be- 
tween the overlap fluctuation, i.e. 60^60iy = 60^ 6^^, we get the following 
correlation 

C(E) = (5) 

where AE is the averaged level spacing. This relation shows that the width 
of the autocorrelation function is directly related to the typical width of the 
individual states. 

In flgure (^, two cases are displayed: in the left part, we illustrate a 
case with F -C AE, which leads to the fragmented response (in some cases, 
this is comparable to the Landau spreading), whereas in the right part, we 
show a typical Ericson case with T ^ AE. On top of this flgure, we show 
the strength function^. As discussed in section 3, instead of performing the 
autocorrelation technique on the strength function, it appears more conve- 
nient to perform it on the flrst derivative of the strength function. Indeed, 
the autocorrelation on the strength function is particularly suitable when it 
is applied to a rather flat average spectrum. However, in general because of 
the global variation of the strength, an average strength has to be removed 
from the distribution in order to extract the fluctuations [0. In the case of a 
resonant average response, we have observed that the width extracted from 

^It should be noticed that, since the energy units can be re-scaled, we have chosen 
arbitrary units to display the various strength functions. 
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the autocorrelation procedure is sensitive to the method used to define this 
average (fitting procedure, averaging interval). This might introduce some 
spurious fiuctuations. A way to remove it, is to consider the derivative of the 
studied spectrum. The derivative of the strength function emphasizes its fine 
structure and avoid the ambiguities on the averaged strength distribution. 
In the middle of figure (0), we show the derivative of the strength function 
and the associated autocorrelation in the bottom. In both cases the width 
of the autocorrelation is directly proportional to the width of the individual 
states. Although the two considered strength function are different, we can 
observe in figure (|1]) that the associated autocorrelation functions are almost 
identical. This comes directly from the fact that the autocorrelation applied 
to the derivative is only sensitive to the smallest scale in the spectra and ig- 
nores long-range correlation. In this figure, since we have imposed the same 
small fiuctuations scale in both examples, the two autocorrelation happen to 
be very similar showing that this technique could be applied on a fragmented 
strength or a Ericson fiuctuation case with the same success. 



2.2 Damping of Collective States 

In order to get a deeper insight into the fiuctuation pattern in the collective 
response and into the information which can be extracted from their study, 
we can now consider the microscopic mechanism responsible for this damping. 
The coupling of the collective motion to the compound nucleus states leads 
to an internal mixing while the relaxation to the continuum is a true external 
decay channel. In reference [|^, the effect of the continuum was carefully 



investigated but no particular assumption on the internal degrees of freedom 
is made. In the present paper we will mainly focus on the complexity of the 
internal mixing. In particular, we will introduce a hierarchy of degrees of 
freedom and couplings. 



2.2.1 Modelisation of the Doorway Mechanism 

Let us first introduce a doorway state \Coll) associated with an unperturbed 
energy Hq \Coll) = Ecoii \CoU) . This state is in general not an eigenstate of 
the total hamiltonian but is coupled to more complexe states. Very often it is 
possible to introduce a hierarchy of complexity between those states. In par- 
ticular, giant resonances are assumed to decay toward two-particle two-hole 
states which are themself coupled to three-particle three-hole states which 
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are damped by more complex states. Such a mechanism can be described 
by considering a basis which can be formed into a hierarchy of states \CoU) , 
ilh)} , {1^2)} , ... at unperturbed energies Hq \in) = Ei^ \in) ■ These states 
allows the definition of more and more complex Hilbert spaces So G Si G S2 
C. . . where Sq contains only \CoU) while Sn includes all the states up to the 
j^th igyg[^ 1^^^ _ gg^]^ Ylow introduce a hierarchy of residual interactions Vi, 
V2, ... where the interaction couples the states of the space Sn^i to the 
states {|in)} • Then, the diagonalization of the Hamiltonian 

Hn = Ho + Vl + ... + Vn (6) 

produces eigenstates |yU„) 

\f^n)=c;^,Coa\Coll)+j:c;„,\^) (7) 

i 

with an energy E^^^ . In the following we will truncate the hierarchy at various 
level and discuss the fluctuations properties of the associated strength. 



2.2.2 The escape width as a coherent decay of the states /i 



At a given level of complexity, the states /i have an infinite lifetime. However, 

in nuclei, eigenstates are coupled to the continuum leading to a finite lifetime 

for each state. This decay to the continuum is also linked to the escape 

width of the collective mode. Indeed, introducing the continuum states 

\k) and a coupling with the states /i through a residual interaction AV^, we 

can estimate the width which is related to the transition matrix elements 
2 

{k\ AV l/i) from the states to the continuum states \k). Using the wave 
'unction , we get 

2 



{k\AV\fi) 



c^^coii {k\ AV \Coll) + c^,, {k\ AV |2) 



(8) 



Since c^^cou and the various c^^j are not correlated, the previous expression 



can be evaluated as an incoherent sum 



{k\AV\^^) 



\Cii,CoU\ 



AV \Coll) 



+ 



E 



I 



{k\ AV \i) , showing that, as far as the decay to the continuum is 



concerned, taking advantage of the Fermi golden rule, we can write 



\Cfi,Coll\ 



(9) 
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Here we have introduced the direct couphng of the collective states to the 

- 2 

continuum: the escape width F' oc {k\ /SV \Coll) . T ^ can thus be split- 
ted into two parts: a collective one r^°" = |c^,co«|^rT = 0^T\ which is 
distributed according to a lorentzian shape (cf eq. ([TI|) ), and an individual 
coupling to the continuum Fj^ = |c^,i|^rl which is not expected to present 
any particular structure. Since the |c^,co//|^ are very small, the might be 
much smaller than F^^. This is actually the case since typical values are the 
eV for F^ and the hundreds of keV for F^ Note that, if we take \c^^coii'^ 
oc where is the number of states coupled to \Coll): N ~ Tcoii/^E 
with AE the average spacing between two states /x, we see that the collective 
contribution to F^ is proportional to F^/A^ = AE /Tcoii which is smaller 
than F^. 

Then one may ask how it is possible to get a large direct decay probability 
for the collective state out of small individual decay rates F^? This is due 
to a coherence effect. Indeed, if we excite the collective state \Coll) and we 
want to compute its decay probability, in terms of the individual probability 
F^, using I Co//) = ccoii,/, 1/^), we get 



F^ 



oc 



{k\ AV \Coll) 



(10) 



This coherent sum can be much bigger than the incoherent analogue to 
J2n \ccoii,fif T/x- This means that it is possible to have a strong direct cou- 
pling from the collective state to the continuum built out of a coherent sum 
of many small decay probability. 

In the following, the direct decay of the collective state is included in the 
individual width F^ of each eigenstates |/i). To do so, the coherence should 
be kept and the width should fulfill equation (^. However, we will often 
neglect the energy dependence of the width F^ in the analytical expressions. 
We have tested numerically that this energy dependence does not affect the 
presented conclusions. 

2.2.3 Microscopic Description of Ericson Fluctuation in Resonant 
Phenomena 



Let us first recall the standard description [11^ of a doorway state which 



includes only the first level. In such a case, assuming (Co//| Vi 
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and a regular spectrum of states with a level spacing Ei^—Ei^_i = AEi, 
the overlap matrix becomes (see appendix 2) 



(vl + vfnyAE!) + (E,, - Ecouf ^^^^ 

which is the standard Lorentzian shape with a typical width V^^^ = 4 
In the limit of a continuous spectrum and of AEi assuming vf/AEi = 
cte, one gets Tcoii = 2nvl/AEi which is equivalent to the standard Fermi- 
golden rule. 

In this derivation no fluctuation in the strength function is introduced. 
This can be traced back to the simplifying assumptions about the density of 
states \i) and the constant interaction matrix elements {CoU\ Vi Let us, 
for example, introduce a fluctuating part of the residual interaction matrix 

2 

element defined by {Coll\ Vi + 6Vi = (1 + Svi-^^). Then the overlap 
matrix can also be splitted into an averaged and a fluctuating part O^^ = 
O^^ + 50^^ , the averaged being identical to the constant interaction case 
(see equation ( [Tl| ) and Appendix 3). Introducing diagonal correlation of the 
fluctuating matrix element 5vi^5vj^ = Cj^^j^ji, we get 

50,^50., ^ 5^,., 501 (12) 

where the expression of S0^_^ is given in appendix 3. Therefore, one get back 
the conditions needed to observe Ericson fluctuations. 

We have investigated this effect using a numerical diagonalization of a 
non- fluctuating Hamiltonian Hi (top part of fig. (^) and of a fluctuating 
one (bottom part of fig. (0)). To describe the finite lifetime of the states 
an imaginary part —i has been added to each energy E"^^. On figure @, 

we can see that a fluctuating residual interaction produces Ericson fluctua- 
tions. The above discussion applies to the second RPA^, |2^, |24 



description of collective excitations which takes into account the decay of col- 
lective states excited through one body operators, i.e. the coupling of states 
built from particle-hole (p-h) type of excitations, into more complex config- 
urations containing two-particle two-hole (2p-2h) states. In this framework, 
one would expect to observe fluctuations of the collective strength related to 
the characteristic lifetime of the (2p-2h) states. 

2.2.4 Fluctuations within a Microscopic description of the lifetime 
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of the decay channels 

In order to get a deeper insight into the fluctuation mechanism, one can 
introduce the second level of complexity 1*2) • The corresponding collective 
response can be computed by introducing as follows: 

• flrst the state \Coll) is coupled to the by the residual interaction 
Vi + 6Vi leading to eigenstates \fii) 

If^i) = cl,Cou\CoU)+Y.cl^,^\^^) (13) 



Then, the states are coupled through the residual interaction V2 
with the states 1^2) in order to built more complex eigenstates |/i2) 

1/^2) = cJ^,,J/ii)+E^i,^J^2) (14) 



Therefore we get 

1/^2) = cl^,^,,cluCoU \Coll) + cj^_,, Ki) + E K2) (15) 

If we assume that the states are independently coupled to different en- 
sembles of states 1^2) through a constant interaction ^2, then, each will 
act as a doorway state toward its own decay channels and so will be spread 
over a typical width F^^ = 4 (f | + f|7r^/Aii^|.) where Aii^2 is nothing but the 
level spacing at the states \i2) coupled to . With this formalism it is 
clear that the overlap matrix presents in fact two scales 

O,, = Olfil^^^ (16) 

the flrst one O^, associated with the flrst level of complexity yUi, and a second 
one 0^2,^1 associated with the spreading of /xi over ^2- Within the presented 
approximation the strength function reads 



_ (17) 

where we have introduced a width V^^ for the states /i2. Describing O^^ as 
in eq. (0) and the 50^^ as fluctuation, we recover to the case of Ericson 
fluctuations as discussed in the previous chapter. 
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2.3 Multiscale Fluctuations and the Doorway Hierar- 
chy 

In the previous section, we have assumed that the doorway state was directly 
coupled to a large ensemble of states with an average intrinsic width. In such 
a case, fluctuations with a typical scale, related to this width, are observed 
in the response spectrum. We have also considered the possibility that the 
damping might be due to a coupling to more complex states. However, we 
have not yet investigated the possibility that this second level of complexity 
might also be fluctuating. 



2.3.1 Multiscale Ericson Fluctuation 

We expect in general fluctuations at this second step of the decay, then the 
strength reads 

^ (^) = ^ E K + 501) E [01.. + ^01..) i^E-E^^Kvi /4 

'jl8) 

Using a lorentzian shape both for O^^ with a width TcoU and for O^^ /^i 
with a width F^^, we can see that two scales F^^ and V^^ are present in the 
fluctuations of the strength function on top of the spreading of the collective 
state Tcoii- This is illustrated in figure |^ (top-right) where we have plotted 
a strength computed with two scales for Ericson fluctuations. 



2.3.2 Fluctuations of a fragmented strength 

Another interesting situation is the case of Landau spreading. Indeed, at 
the RPA level of description P^, the strength often appears splitted in 
several components |Co//„). In such a case, at least three scales might be 
identified. The Landau spreading can be directly observe as the fragmenta- 
tion of the strength while the width of each fragment can be identified with 
the lifetime of each individual collective state. Moreover, as discussed above 
the strength presents fluctuations on top of each resonant line-shape which 
are characteristic of the various decay channels widths. Indeed, if we assume 
that the different fragments are decaying toward independent complex states, 
the strength is simply the sum of the various components associated with a 
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collective excitation |Co//„) 

s{E) = y: I {Coll \coiin)f ( E -ll^}£^Ml\ (19) 

Such a fragmented strength is illustrated in figure || (top-right). 



3 Observation of fluctuations at different scales 

The observation of the collective state width Tcoiu the largest scale in the 
strength, is rather straightforward using standard techniques such as fits, 
variance estimation or even autocorrelation calculations. The extraction of 
characteristic widths from a spectrum presenting many different scales of 
fluctuations is a more complicated problem. In this chapter, we discuss sev- 
eral methods and propose new approaches for extracting signals of multiscale 
fluctuations. 



3.1 Standard method and its extension. 

We have seen in the previous chapter that the autocorrelation function is a 
useful tool when fluctuations and fragmentation of the response are involved 
(see flgure (|1|)). On the one hand, when the autocorrelation techniques is 
used directly on the strength distribution, the obtained signal can be related 
the total width of the collective mode. A standard technique to extract the 
fluctuation properties of a flne structure on top of some smooth "background" 
is to subtract an averaged distribution |]TB[ from the signal. However, the 



obtained result appears to be dependent on the method and parameters (such 
as the smoothing interval) used to deflne this average strength. In order to 
overcome this problem, we have performed the autocorrelation analysis on 
the derivative of the strength distribution. The derivative is very sensitive 
to the fluctuations over the smaller scale present in the spectrum since they 
are associated with very rapid variations. In this case, the autocorrelation 
function gives half of the width of the smallest fluctuations scale. 

These two types of autocorrelation analyses, on the strength and on its 
derivative, are illustrated in flgures (131) for two multiscale (two scales) 

Ericson fluctuation (left of flg. (|^)) and a fragmented strength with flne 
structure (right of flg. (|)). 
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In both cases (multiscale Ericson of fragmented), the largest and the 
smallest scales in the fluctuations can be extracted. Considering now the 
fragmented strength case, an intermediate scale can also be observed. How- 
ever, when Ericson fluctuations occurs at many different scales, the auto- 
correlations techniques applied to the strength and its derivative only gives 
access to the largest and the smallest scales of the fluctuations and interme- 
diates scales could not be obtained. 

During the past years, different alternative techniques have been proposed 
to extract properties of fluctuations in the decay of collective motions ||2^, 
pO|] . In particular, in a second RPA picture, where only two-body 
correlations are retained, a possible fractal nature of the fluctuations has 
been discussed [^, When higher correlations are considered, we 

have shown that a hierarchy of well-separated scales may also be present. In 
the next section, we will propose a novel technique, that may help to get 
either the well separated scales or the fractal nature of nuclear decay. 



27, 
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3.2 Entropy index for multiscale fluctuation 

Recently, a method based on the definition of an entropy index has been 
proposed by Hwa[^] in order to extract scaling behavior in fluctuating signal. 



The entropy index appears as a good indicator of the existence of different 
scales in a strongly fluctuating spectrum. We have adapted this method to 
the nuclear response case. 

The total energy interval /S.E = Em ax — Emin, is first divided into n 
bins of resolution 6E {n = AE/SE). The signal in each bin is analyzed 
according to a simple overlap with a step function changing from —1 to +1 
in the middle of the considered binQ. Therefore, at every scale 6E, in each 
of the bin j, a coefficient Dj {6E) is defined as 



Dj {6E) = dE S {E) sign {E - (j - 1/2) 6E) (20) 



The coefficients Dj can be considered as a coarse grained derivative of S. In 
order to focus on the global properties of the fluctuations at a given size 5E 
, an entropy factor K (SE) can be defined as 

K {5E) = -Y.^i (^^) log Wj (SE) (21) 

^ j=l,n 



*This method can be viewed as a wavelet analysis. 
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where the coefficient Wj (SE) = Dj/ (Dj) are nothing but the coefficients Dj 
normahzed to their averaged value ((-Dj) = l/nJ2j=i,nDj)- 

In the reference ||31|| , it is shown that a hnear decrease of K [6E) as a 
function of log (SE) characterizes the existence of fluctuations at all scales. 
When fluctuations with specific scales are considered, a different behavior is 
expected. Having in mind that the coefficients Wj [SE) might be interpreted 
as a normalized coarse-grained derivative, K {6E) should remain almost con- 
stant between two typical scales since in-between these scales, the coarse 
grained derivative of S do not vary much. However, when going from one 
scale to another, the derivative looses part of its structure and one is ex- 
pecting an entropy variation. This situation is illustrated in figure @ for 
three different cases, where respectively fluctuations over one, two and three 
scales have been introduced in the strength. Note that, in order to avoid 
problems due to the limited number of bins for large bin size, we have con- 
sidered a function defined as a repetition of the strength function instead 
of the strength itself. In practice, we used 31 repetitions of the strength. 
The entropy index is then applied on the new function inside a large energy 
interval, enabling to have a large number of bins in the energy region under 
interest {6E < AE). From fig. (^, one can observe that the evolution of 
K (SE) indicates the presence of respectively one, two and three scales by a 
change in its curvature. In order to emphasize this evolution, we have plot- 
ted in figure (|^) a numerical estimation of the second derivative K" (SE) of 
the entropy index. The presence of one, two or three scales is thus signed 
by the presence of respectively one, two and three minima in these second 
derivative. 

In all the presented cases, K{6E) is a good indicator of the different char- 
acteristic scales in the nuclear response. Using the positions of the curvature 
variations, the various scales can be roughly estimated. 

In this section, the entropy index method was applied to model where 
scales are well separated. This situation may append at rather low excita- 



tion energy where few degrees of freedom are coupled to collective states|32 



When excitation energy increases, the number of internal and external de- 
grees of freedom involved in the coupling becomes very large [ 33 1 and a sta- 
tistical treatment is required []T^ |TB|, In this case, many different scales 
are expected and can be uncovered through the entropy index method. This 
method can also be used in more complex situations where mixing of the dif- 
ferent scales is expected. It can even sign self-similar fluctuations |^, |2^, |29 |. 
It should be noticed that the smallest scales expected in the strength are re- 
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lated to the lifetime of the compound nucleus states which are less than the 
eV. However, experimental resolutions have not yet reached this degree of ac- 
curacy |3^ . Elaborated techniques based on statistical assumptions have 
been applied in order to extract the properties of invisible fine structures 



The entropy index method can be viewed as a model independent way to 
extract information about fluctuations. Only scales larger than the exper- 
imental resolution can be accessed but no assumptions on the statistical 
properties of the studied spectrum are needed. 



4 Critical discussion of the Doorway hierar- 
chy picture 

In order to observe multiscale fluctuations, we have assumed a particular hi- 
erarchy in the Hamiltonian following the general scenario of a gradual com- 
plexification of a collective motion until it reaches the compound nucleus 
chaotic states. It is in fact the implicit assumption of many simulations such 
as the extended mean-field approach which takes into account the two-body 



collisions as a damping mechanism p^, |37| , p8| , p9| , |40[| . In these hier- 
archy of couplings is implicitly assumed, the main ansatz being that cutting 
this hierarchy at any level leads to an approximate strength which is only 
slightly modified when the next level is introduced. From a quantal point of 
view, this means that an ensemble of states ordered by increasing complexity, 
such as states with increasing number of quasi-particle excitations, can be 
defined and that each level of complexity can be considered as a perturbation 
on top of the previous level. However, there is no a priori reason that such a 
robust hierarchy exists. Indeed, in quantum mechanics, a modification of few 
matrix elements is often sufficient to introduce an important rearrangement 
of the whole spectrum. Therefore, even if the hierarchy of doorway seems 
valid for several decay steps, it is always possible that the introduction of the 
next level of complexity deeply transforms the overall picture. For example, 
the coupling of overlapping resonances through their common decay chan- 
nels has been already discussed in connection with nuclear relaxation (see 
for instance [0, |15[). Specific coherent effects, like the Dicke supperradiance 
well known in quantum optic^], could also be present in nuclear spectra 



and could suppress Ericson fluctuations|]T3[. Such coherence effects can also 
affect the nuclear response when a hierarchy of intrinsic degrees of freedom 
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is considered. 

Let us consider a single collective doorway state \Coll) damped through a 
residual interaction Vi and V2 assuming that {Coll\ Vi = vi, \ V2 \i2) = 
V2 and that the second level of complexity 1^2) is characterized by a level 
spacing AE2. Then we can take the limit f | — > and AE2 keeping 
Fi = 27rf|/A£'2 constant. In the case of a regular ensemble of states ii with 
a spacing AEi, we can write the overlap matrix (see Appendix 4) 

o =|pM2 |2 = 

^-^^2 \^Coll\ p2 r, AE? , / P \ -n \2 

^ {E^, - Ecouf + ^ {{E,, - Ecou) tan (vrl^) - ^) 

(22) 

where we have introduced Tcoii = ^nvf/AEi. When Ti ^ AEi. In such a 
case, the overlap has a typical scale of fluctuation not related to Fi but 
to AEi since the denominator goes to infinity each E^^ = (n + l/2)A£'i in 
eq. (0). 

Again we have tested that all type of fluctuations of the residual Hamilto- 
nian as illustrated in figure (^. This figure presents strong oscillations with 
a width related to AEi as shown by the autocorrelation function. Therefore, 
the occurrence of interference like pattern characteristic of the level spacing 
AEi should be considered as generic. We would also like to mention that a 
strong coupling between the first and second decay channel induces a nar- 
rowing of the apparent width of the collective states. This effect is analogous 
to the motional narrowing described in[^|. 



This illustration stresses the fact that in the case of a coupling of several 
states through their decay channels one should expect interference pattern 
due to a feedback of the decay channels on the properties and the coupling 
of the doorway states. Then, the fluctuation pattern can reflect the typical 
level spacing and not the typical width. 



5 Conclusion 

In conclusion, we have studied the role of several steps in the decay of door- 
way states We have shown that, in some cases, in particular when the decay 
channels of each states encountered at each level of the decay cascade can 
be considered as independent, one should expect to observe several scales 
in the fluctuation pattern of the response function. In such the var- 

ious fluctuations are characteristic of the typical lifetime at each level of 
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complexity encountered in the path from the doorway toward the chaotic 
compound nucleus states. We have shown that this multiscale Ericson fluc- 
tuations can be partially revealed by applying the autocorrelation technique 
on the strength. However, standard methods which involve diff^erent smooth- 
ing procedures seem to bias the obtained results. In order to overcome this 
problem, we have proposed two new techniques: the first one consists in 
computing the autocorrelation function on the derivative of the strength dis- 
tribution in order to extract the smallest scale in the fluctuations, the second 
one, called Entropy index, gives information on how many scales does exist 
in the fluctuation pattern. Therefore, high resolution experiments are indeed 
very interesting tools to study the damping mechanism of collective states. 

Finally, the analytical results as well as the numerical simulations of the 
last chapter illustrate that a hierarchy of scales in the fluctuation, directly 
related to the lifetime of each level of complexity, is not the generic case 
and that often the fluctuations can be related to other physical quantities. 
In particular we have investigated the effects of decay channels interacting 
via the states encountered at a higher level of complexity. In this case, 
interferences between different decay channels are found and the induced 
fluctuations are not related to the decay width of the considered states but 
to their spacing. In such a case, experiment with a very high resolution will 
present fluctuations which are characteristic of the level spacing and not of 
the width as in Ericson fluctuations. 

Caution should thus be used when interpreting the observed width of 
the fluctuations. Indeed, the interferences, due to an interaction of some 
states sharing the same decay channels, may produce narrow fluctuations 
which may mimic long lived systems. This is an important flnding which can 
modify the interpretation of Experimental data. 
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APPENDIX 1: Ericson Fluctuations in doorway processes 

Let us consider the spectral response 
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When we introduce an averaged overlap and a fluctuating one, | (/i \Coll) |^ = 
O (E^) + 60^, the response can also be split into two components. On the 
one hand, it contains an averaged response function R (E) = J^^ dE' R {E') 
which can be written as for E > {] 

R(E)^^ f dE' p (E') -pOiE) (24) 

where we have introduced the density of states p{E) = J2ij,^ {E — E^), the 
average density of states p = J^^dE' p{E') /SE = 1/AE and an averaged 
width r (the integration interval 6E being taken much larger than F). On 
the other hand remains a fluctuating term 

These fluctuations still contain information about the averaged width T as 
it was demonstrated by Ericson|T^. Indeed, the autocorrelation function of 
this noise is given by 

1 rEo+SE 

Cr{E) = — dE' 5R* {E') SR {E' + E) (26) 

oE JEo 

Cn{E) = y— dE'—^ (27) 

' jn^dEJEo E'-S; E' + E-S, ^ ^ 

where we have introduced complex energy notations £^ = E^ — If we 

assume that 6E ^ F, we can extend the integration to inflnity, so that, taking 
advantage of the relation (x — y)~^ [x — z)~^ = {y — z)~^ {^{x — y)~^ — {x — zY 
it is possible to perform a Cauchy integral leading to 

If we assume random correlations between the overlap fluctuation 

80^80, = 801 V (29) 
in average, we get the correlation 



47r 80"^ 

Cr{E) = , (30) 

^ ^ AE E^ + T 
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A similar equation holds for the strength function associated to equation 

Indeed, we can define the autocorrelation function of the strength distribution 

as 

1 r Eo+SE 

C{E) = — dE' 5S {E') 5S {E' + E) (32) 

oE JEo 

Noting that S{E) = -l/(2i7r) {R{E) - R*{E)) we have 



C{E) = I ^ . (33) 

This relation shows that, even in the case of strongly overlapping states, one 
can still extract information about the typical width by looking at the fluctu- 
ating part of the considered autocorrelation spectrum as it was demonstrated 



m 



APPENDIX 2: Simple Doorway Picture of a Damping Mecha- 
nism 

The diagonalization of the Hamiltonian Hi = Hq + Vi produces eigen- 
states l/ii) associated with the eigenenergies E^^ which fulfill the following 
dispersion relation 

E^,-Ecou = h{E^,)vl (34) 



where we have considered a constant interaction vi = {Coll\ V 
overlap matrix O^^ = \ {fii |CoZ/)|^ reads 



. The 



^'^^ = 1 , f rV (35) 
J- + 72 [E^JVi 

In the previous equations, we have introduced the functions 

n [E - Ei^) 
which are related by a recurrence relation 

fn+i [E) = ~n (37) 
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In order to get the usual Lorentzian shape, we assume a regular spectrum 
of energy Ei-^: Ei^ — -Ej^-i = A£'i with ii running over both positive and 
negative integer values. Then, using the relation I]j 1/ ~ = ^ cot (vrx) , 
the dispersion relation becomes 

- £co„ = /■ ) = ^ cot (Sfl^^MJ (38) 

On the other hand, using the relation ( |37D and taking advantage of the 
dispersion relation (|38|) , /2 (-E/^J can be recast as 



/.(g„)-^^+ '^'";f"°"^' i (39) 

so that the overlap matrix simply reads 



?;2 

O = -i (40) 

{vl + vW/^El) + {E^,-Ecoiif ^ ' 

which is the standard Lorentzian shape with a typical width Vjj^^ = 4 (t;^ + v\tt'^ / AEf). 

APPENDIX 3: Microscopic Description of Ericson Fluctuation 
in Resonant Phenomena 

Following Appendix 2, let us, for example, introduce a small fluctuating 
matrix 6Vi. Then we define a fluctuating part for the /„ functions 

^ ^ 2 

where we have used the notation {Coll\ Vi + 6Vi = + SviJ . The 
new dispersion relation reads 

E,, - Ecoii = vf{h iE,J + 5f, (E,,)) (42) 

while the overlap becomes 

= l+WHf2{E,,)+Sh{E,J) ^^^^ 

In the case of a regular spectrum of states , using the relations derived 
above, the dispersion relation becomes 

E,, - Ecou = vf cot (vr^^^^i^;^) + 5/i {E,^)^ (44) 
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so that, we may use this relation in order to compute /2 (-E^i) 

/2 (E,,) = + v-^-' (e,, - Ecoii - vf5h {E,Sf) (45) 

The latter equation can be introduced in the overlap matrix leading the 
perturbative result 

O^, = O^, + 50,, (46) 

with an averaged value, O^j, identical to the one derived in appendix 2 and 
with a fluctuating part given by 

50,, = 01, (2 {E,, - Ecoii) 5h (E,,) - vf5h (E,,)) (47) 
We finally get : 



<^O^i=Oj^EKHi(^n) (48) 

with 



Ml 



So (E)^ - r^jE^^-Ecoii) 1 \ 

This is a fluctuating correction to the averaged overlap which looks hke the 
one needed in order to observe Ericson-like fluctuations. Indeed, the correla- 
tion 50,,50u, reads as a function of the correlation of the fluctuating matrix 
element 5vi,5vj, — Ci,j, 

50,,50., = 01,01, Y: CnnSo,, (E,,) 5o,, (Ej,) (50) 
hji 

If we now assume uncorrelated fluctuations for the matrix elements Ci,j, — 
Ci,5i,j, we get 

50„50., = Ol,0l,Y^c„5o„ {E,,) 5o,, {E,,) (51) 

because of the dispersion relation, the energy E„ is always very close from 
one particular state which we note associated to an unperturbed 
energy E^n . Then, in 50,, , mainly this state i^^ is participating to the sum 

over ii so that 50„ 0'^„^v^no„ {p^n^ and the correlation reads 

50„50., « 5,,., 01, c,.i 5ol^ (e,.. ) = 5,,^, 501 (^2) 
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Therefore, one get back the correlations needed to observe Ericson fluctua- 
tions. However, it should be noticed that compared with the simple constant 
correlation assumed in eq. (^) one may have here a smooth energy de- 
pendence of SOf^^. This however do affect the conclusions reached about the 
width of the autocorrelation function. 

APPENDIX 4: Interferences in the Damping of a Resonance 

We can finally consider the related case of a single doorway collective 
state \Coll) damped through a residual interaction Vi toward a first level of 
complexity described by the states which are then coupled to many states 
1^2) by a second part of the residual interaction V2. Then, the diagonalization 
oiH = HQ+V1+V2 (where \Coll), and |z2) are eigenstates of f^o) produces 
eigenstates 

(53) 

The Schrodinger equation leads to 

h 

Cf,,,n (^M2 -Ei,) = C^^ coll (^ll \Coll) + J2 ^f.2,i2 (^ll '^2 K2) = VlC^,,CoU + ^ 



S2,i2 (^M2 - Ei^) = ES2,n ^2 = V2al 



«2 

), 

M2 



where we have assumed that {Coll\ Vi = Vi, \ V2 K2) = "^2 and where 
we have introduced a^^ = c^2,n' ^^2 = ^^2 '^/i2,i2- '^^^ ^^^^ equation gives 



''li2,CoU 771 771 

^^l2 ~ ^Coii 
so that the two other relations reads 



^n2,Coii — TP (57) 



M2,2l 



M2,«2 



/^2 

V2a 



Ei^ \ E, 



M2 



E, 



+ V2a 



Coll 



M2 



^2 



-^M2 - -^12 



(58) 
(59) 
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Therefore, the following relations between Oi and 02 can be expressed as 

< = E /'!^ ^v,alf,{E,,) (60) 



/ ^ TTi TP TP TP TP TP \ TP TP V'Vf'S/ 



where we have introduced the functions /„ and F„ 

= ^^-1^ (62) 

/„(E,J = Y. fj, \ (63) 

Then the eigenenergy En^ fulfills the following dispersion relation 



2 



1 = J, + ^2 /i (^.2) ) F, (E,,) (64) 

\ J^ti2 ~ ^CoU ) 

Moreover, the a coefficients can be found using the normalization condition 

21^ /21 \2 



— ,p p ,2+^2<2/2(-g.2)+ U F ^2(^/^2) 

(E^2 - Ecoll) V^M2 - ^Coii / 

(65) 

which leads to a condition on a, through the use of the relation (§^), in order 
to remove a^^? 

1 = ,^ ""'t' ,2 +^2 <2 /2 (^.2)+ f + -2^ < /l (^M2)) ^2 (i?,2; 

(E^2 - Ecou) V^M2 - ^CoZZ ) 

(66) 

allowing to extract the coefficient a\ 

<2 = f ,^ ,2 + ^2 /2 (i?M2) + + ^2 /l (^.2) V ^2 {E,.,] 



{Ef,2 - Ecou) \Ef,2 - Ecoii 



2 

Then the overlap matrix, = |c^2,CoZ/| , reads 



(67) 



(68) 
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Let us now assume as usual that the second level of complexity \i2) is char- 
acterized by a level spacing A£'2, then we can write for fi 

ME) = -^cot(n-^] (69) 



AE2 \ AE^_ 

Using the relation ([37|), the corresponding relation for /2 can be obtained 



We get the dispersion relation 

then we get for the coupling coefficient cf^^ 
2 

O 



M2 



t>of 



2 2 
2^1 



\ * i v"''2 y / 

(72) 

where we introduce Fi = ^nv^j !\E2- If we assume that the interaction 
f I — i> we can approximate the overlap by 

Q _ |„/^2 |2 _ /^^N 

^1 /'P P ^2 _L ( ^n-^Coll „,2 



Let us now consider the case of a regular ensemble of states ii with a spacing 
AEi, we can write for Fi 



F.(B„J = ^cot(.^) (74) 



and the overlap reads 
O 



'M2 — \^Coll\ 



,2 



^ (E^, - Eco.)' + {{E,, - i?co.) tan (vrif ) - ^ 

(75) 

where we have introduced TcoU = 



24 



References 

[1] G. F. Bertsch, P.F. Bortignon and R. A. Broglia, Rev. Mod. Phys. 55 
(1983) 287. 

[2] G.F. Bertsch and R.A. Broglia, Oscillation in finite quantum systems, 
Cambridge (1994). 

[3] Ph. Chomaz and N. Frascaria, Phys Rep. 252 (1995) 275. 

[4] P. A. Moldauer, Phys. Rev. 171 (1968) 1164 

[5] P. Kleinwacher and I. Rotter, Phys. Rev. C32 (1985) 1742. 

[6] S. Drozdz, S. Nishizaki, J. Speth and J. Wambach, Phys. Rep. 197 
(1990) 1. 

[7] Ph. Chomaz, NucL Phys. A569 (1994) 203c. 
[8] Ph. Chomaz, Phys. Lett. B347 (1995) 1. 

[9] B. Lauritzen, P.F. Bortignon, R.A. Brogha and V.G. Zelevinsky, Phys. 
Rev. Lett. 74 (1995) 5190. 

[10] V. Zelevinsky, Annu. Rev. Nucl. Part. Sci. 46 (1996) 237. 

[11] J. Winchenbach, K. Pingel, G. Holzmarth, G. Kiihner and A. Richter, 
Nucl. Phys. A410 (1983) 237. 

[12] T. Ericson, Phys. Rev. Lett. 5 (1960) 430; Ann. Phys. 23 (1963) 390. 

[13] C. Mahaux and H.A. Weidenmiiller, Shell Model Approach to the Nuclear 
Reactions, North Holland, Amsterdam (1969). 

[14] V.V. Sokolov and V.G. Zelevinsky, Nucl. Phys. A504 (1989) 562; Ann. 
Phys. 216 (1992) 323. 

[15] V.V. Sokolov, I. Rotter, D.V. Savin and M. Miiller, Phys. Rev. C56 
(1997) 1031; Phys. Rev. C56 (1997) 1044. 

[16] P. Ring and P. Schuck, The nuclear many-body problem, Spring- Verlag, 
New- York (1980). 



25 



[17; 

[18, 

[19 
[20 
[21 
[22 
[23 

[24 
[25 

[26 

[27; 

[28; 

[29 
[30 

[31 
[32 
[33 



A. Bohr and B. R. Mottelson, Nuclear Structure, vol. I (Chap. 2). 

G. Kilgus, G. Kiihner, S. Miiller, A. Richter and W. Kniipfer, Z. Phys. 
A326 (1987) 41. 

V.V. Sokolov and V.G. Zelevinsky, Phys. Rev. C56 (1997) 311. 

H. Suhl and N. R. Werthamer, Phys. Rev. 122 (1961) 359. 
J. Sawicki, Phys. Rev 126 (1962) 2231. 

J. Da Providencia, Nucl. Phys. 61 (1965) 87. 

C. Yannouleas et al, Nucl. Phys. A397 (1983) 239. 
C. Yannouleas, Phys. Rev. C35 (1987) 1159. 

C. Yannouleas and S. Yang, Nucl. Phys. A455 (1986) 40. 

D. J. Rowe, Nuclear Collective Motion: Models and Theory , (Methuen, 
London, 1970). 

S. Drozdz, S. Nishizaki and J. Wambach, Phys. Rev. Lett. 72 (1994) 
2839. 

S. Drozdz, S. Nishizaki and J. Wambach and J. Speth, Phys. Rev. Lett. 
74 (1995) 1075. 

S. Drozdz, A. Trellakis and J. Wambach, Phys. Rev. Lett. 76 (1996) 
4891. 

A. Bershadskii, Phys. Rev. C59 (1999) 3469. 



H. Aiba and A. Matsuo and , Preprint [nucl-th/9904037| , to be published 
in Phys. Rev. C. 



R.C. Hwa, [physics/ 9809"04T 



D. Lacroix, Ph. Chomaz and S. Ayik, Phys. Rev. C58 (1998) 2154. 

N. Frazier, B.A. Brown and V.G. Zelevinsky, Phys. Rev. C54 (1996) 
1665. 



[34] V.V. Sokolov and V.G. Zelevinsky, Phys. Lett. B202 (1988) 10. 

26 



[35] S. Kamerdzhicv, J. Lisantti, P. von Neumann- Cosel, A. Richter, G. Ter- 
tychny and J. Wambach, Phys. Rev. C55 (1997) 2101. 

[36] S. Ottini, N. Alamanos, A. Drouart, W. Mittig, H. Savajols, private 
communication. 

[37] C.Y.Wong and H.H.K. Tang, Phys. Rev. Lett. 40 (1978) 1070; Phys. 
Rev. C20 (1979) 1419. 

[38] R. Bahan, Y. Alhassid, and H. Reinhardt, Phys. Rep. 131 (1986) 1. 

[39] Y. Abe, S. Ayik, P. G. Reinhard and E. Suraud, Phys. Rep. 275 (1996) 
49. 

[40] S. Ayik and M. Dworzecka, Phys. Rev. Lett. 54 (1985) 534; NucL Phys. 
A440 (1985) 424. 

[41] R.H. Dicke, Phys. Rev. 93 (1954) 99 

[42] R. Brogha, T. D0ssing, B. Lauritzen and B. R. Mottelson, Phys. rev. 
Lett. 58 (1987) 326. 



27 




-40 -20 20 40 -40-20 20 40 
E (arb. units) E (arb. units) 




1234012345 
E (arb. units) E (arb. units) 



Figure 1: Top: Illustration of the strength computed for a collective excita- 
tion coupled to an ensemble of states (decay channel). In left part, few states 
with a width F = 3 smaller than the spacing {AE ~ 10) is used. In right 
part, the opposite case of a dense spectrum of decay channels {AE ~ 0.5) 
with an average width F = 3 greater than the spacing exhibit typical Eric- 
son fluctuations. In both cases, the collective state is assumed to be spread, 
giving a variance around 30 in the considered interval. Middle: derivative of 
the strength function. Bottom : auto-correlation performed on the deriva- 
tive of the strength function. In both cases, the averaged width of the decay 
channel F can be extracted. Note that half of the width is measured on the 
autocorrelation of the derivative of t8§ strength function. 
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Figure 2: Illustration of the strength computed for a collective doorway 
state coupled to a series of decay states damped with a characteristic width 
larger than the spacing. Top: when the coupling matrix elements are not 
fluctuating. Bottom: when the residual interaction contains random terms. 
In both cases, the total width is roughly equal to 30 whereas the individual 
width and spacing are F^j = 3, AE'i = 0.5. 
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Figure 3: Top: Strength distribution presenting multiscale fluctuation. Left: 
Ericson fluctuations over two scales (equivalent to flg. 4). Right: Small Er- 
icson fluctuations on top of a fragmented strength (equivalent to flg. 5). 
Bottom: associated derivatives of the strength distributions. In both pre- 
sented cases, the total spreading width is equal to Tcoii — 30 whereas the 
scales of fluctuations of the two decay channels are F^j = 3 and F^j = 0.5. In 
each flgures, the associated autocorrelation function respectively performed 
on the strength (top part) or on the derivative of the strength (bottom part) 
are presented in insert. 
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Figure 4: Top: from left to right, strength function with respectively one, two 
and three scales of fluctuations (F = 30, 3 and 0.5) are displayed. Bottom: 
Evolution of K (SE) as a function of the bin size 6E for the three considered 
strength function with one (circle), two (diamond) and three scales (cross). 
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Figure 5: Numerical estimation of the second derivative (noted K" (SE)) 
of the entropy index as a function of the bin size SE. K"{()E) is associated 
with strengths where respectively one (top), two (middle) and three (bottom) 
scales of fluctuations are present. In each figures, the scales of fluctuations 
are indicated by thick vertical lines. 
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Figure 6: Illustration of the strength computed for a doorway state cou- 
pled to many decay states themselves coupled to an ensemble of states and 
computed introducing all types of fluctuations in the residual hamiltonian. 
For the calculation, we have taken Tcoii — 100, AEi = 3 , Fi = 30 and 
AE2 — 0.5. In addition, a small width F^j = 1 has been assumed for the 
eigenstates of the hamiltonian in order to draw the strength function. In 
insert, the autocorrelation function is plotted demonstrating that the widths 
are proportional to AEi. 
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